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Within a phase field approach which takes the strain-induced elasticity into account, the kinetics 
of the coherent order-disorder transition is investigated for the specific case of Al^Zr alloy. It is 
£SJ ■ shown that a microstructure with cubic LI2 precipitates appears as a transient state during the 

decomposition of a homogeneous disordered solid solution into a microstructure with tetragonal 
DO23 precipitates embedded into a disordered matrix. At low enough temperature, favored by a 
• weak internal stress, only LI2 precipitates grow in the transient microstructure preceding nucleation 

of the DO23 precipitates that occurs exclusively at the interface of the solid solution with the LI2 
precipitates. 

Analysis of microstructures at nanoscopic scale shows a characteristic rod shape for the DO23 
precipitates due to the combination of their tetragonal symmetry and their large internal stress. 
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I. INTRODUCTION 



Macroscopic properties of materials depend to a large extend on the microstructures they present at a mesoscopic 
scale. For the case of alloys the most efficient way to control the formation of microstructures is by phase transfor- 
mations. Well known prototypes are the 7 — 7' superalloys, as Ni — Al (see [jlj). The microstructure of such binary 
alloys consists of stable ordered domains, dispersed in a disordered fee matrix. The symmetry of the ordered domains 
is called LI2: the atoms of the minority species are placed at the corner of the fee motif and the atoms of the majority 
, O , , species occupies the center of the faces of the same motif. 

The macroscopic properties are controlled by the size of precipitates, their spatial distribution, and their ability to 
t-h , resist to coarsening. In this context, we present a theoretical study of the dynamics of microstructures in an aluminum 
t> ■ based alloy, namely Al — Zr. 

In this system, for low enough concentration, the equilibrium ordered structure is the tetragonal DO23 phase. Its 
motif is obtained from LI2 structure with anti-phase boundaries in (100) directions with the period of 3. In Al — Zr the 
cubic LI2 phase is known to be metastable at all temperature T . This is confirmed by ab-initio electronic calculations 
at T = OK m and |?J which show that the energy difference between Ll 2 and DO23 is about 0.863 x 10 8 J/m 3 (or 
' 9 AmeV / atom) , in favor of DO^j,- However, the lattice misfit of DO23 with respect to pure Al is significantly larger 
than LI2. Hence, the interplay between the chemical energy and, for coherent microstructure, the elastic energy may 
induce various precipitation processes. 

The aim of the present work is to investigate the decomposition processes and to analyze the resulting microstruc- 
tures at mesoscopic scale. In that aim, we use a phase field approach, where the incoherent chemical energy is 
represented by a Ginzburg-Landau free energy supplemented by a strain-induced elastic energy, in the form proposed 
by Khachaturyan Using numerical simulations, we show that elastic stress favors the LI2 order at low Zr 

concentration and low temperature and in that range the kinetics generates LI2 rather than DO23 precipitates. If 
£j , the coherent elastic interaction is not included, only DO23 precipitates nucleate. 

Our simulations prove that one may obtain metastable LI 2 phase in specific conditions that correspond to the 
. £^ \ experimental conditions leading to a LI2 microstructure. The lifetime of that microstructure may be infinite compared 
to the simulation time, provided the temperature is not too high. Once the LI2 precipitates have grown at sufficiently 
low temperature, the system may be brought to a higher temperature, where no LI2 would nucleate if the system 
had been quenched directly to that temperature. If the already developed LI2 microstructure is aged again at higher 
temperature, then the DO23 precipitation starts and DO23 precipitates nucleate exclusively along the interfaces 
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between LI2 inclusions and the matrix. Once DO23 precipitates are formed, they grow and consume LI2 domains. 
The latest transformation cannot be reversed with decreasing temperature. 

In Sec.|J, the general principles of the kinetic model are presented. In Sec. Ill, we outline how the microelasticity 
contribution can be included into the model. The implementation of the phase transition kinetics, the results obtained 
and their interpretation are given in Sec.|^. Conclusions and perspectives are drawn in Sec.[v|. 



II. GINZBURG LANDAU FUNCTIONAL 



Our mesoscopic method is a time-dependent Ginzburg-Landau kinetics driven by a functional with two parts: 
first a Ginzburg-Landau functional that includes the chemical interactions and the interface effects and second a 
strain-induced elastic energy. 

The explicit form of the Ginzburg-Landau (GL) free energy is imposed by symmetry rules. The first step consists of 
identifying the long range order (LRO) parameters that represent the ordered phases we want to study. In the present 
situation, we should introduce the LRO parameters of DO23 and Ll2- Nevertheless, to simplify the writing of the GL 
free energy, we choose to replace the DO23 LRO parameters with the DO22 ones as these phases are both tetragonal. 
What differs between them is only the periodicity of the anti-phase boundaries that is 2 for DO22 instead of 3 for 
D023- As our aim is not to investigate the competition between the latter phases, our description only requires to 
capture the tetragonality of the possible stable ordered phase. 

The LI2 phase is simple: it consists in three independent parameters r/i, ry 2 , rj 3 which correspond to the amplitudes 
of the three waves that contribute to Ll2- The waves correspond to the three vectors of the reciprocal space K\ — (100) 
, K2 = (010), K3 — (001), respectively. If the microstructure consists only in LI2 domains embedded into a disordered 
fee matrix, the local concentration c(R) would be given by 

c{R) = c {R)+ VjifyexpiiZnKj.R) (2.1) 

i={1.2,3} 

where the quantities co(i?) and r)j(R) vary slowly in space. For a DO22 structure, the probability c(i?) to find a atom 
Zr at position R can be written as follows: 

c(R) = co(R) + exp ( i2nK j- R ) + lj( R ) ex P (iZnQj-R) + 

J'={1,2,3} 

7*(i?)exp(-i27rQ i .i?) (2.2) 

where 7* is the complex conjugate of the amplitude parameter jj and Q\ = (1/2 1 0), Q2 = (0 1/2 1) and Q3 = 
(1 1/2) blong to the reciprocal cubic lattice. This choice is not unique but is sufficient. For example, (1/2 1) can 
be replaced by (1/2 1 0)); each index j correspond to a possible orientation of the tetragonal transformation yielding 
DC>22- For each of the three orientational variants there are four translational variants. For example, for a perfect 
DO22 phase with the orientation j = 1, one has four translational variants defined by either 7* = rji, 7* = —r)i 
7J = —ir/i 71 = ir\\. For simplicity we drop the complex variants. It does not affect the description as there are still 
two translation variants for each orientation of the DO22 phase. 

We note that the state 7^ = and r/j ^ for j = 1,2,3 leads to a density c(R) that describes a LI2 phase. A 
perfect LI2 implies that the three (100)-type waves have the same amplitude. As LI2 preserves the cubic symmetry 
there is no orientational variant but only four translational variants. 

We now develop a uniform Landau functional as a function of {cq, r\ 3 ;, 7^}. The terms of this analytical function 
/i(co, f]-, 7) = Fq./l are selected to fit with the symmetry of the fee lattice, that is each term must be invariant under 
any operation of the fee symmetry group. Formally there is no other rule to realize a polynomial expansion but 
the simplest form is probably the best. In practice, other conditions (see below) must be satisfied. We propose the 
following adimentional functional: 

f L = A n (c - ci) n + ^-{c 2 - c ) Y Vi - Mmrnvd, 

i 

j=i±l 

+f (c 2 ~ c oh f - ^ 4 + f £ 7 4 2 7] + f if (2.3) 

j=i±l 
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where order parameters indices are written modulo 3. The hat symbol points out the adimensional quantities: here 
Jl is an adimensional free energy density. The sets of coefficients {Aj}, {-B?}, {Kj} and {01,02} should be function 
of t he t emperature but are chosen as constant parameters for simplicity. The first term in the right hand side of 
Eq.( [2.3| ) is the disorder contribution as it does not depend on the order parameters. The power n of that term may 
have two values n = n + = 2 and n = n_ = 8 depending whether cq > c\ or not. We introduce that property to adjust 
the topology of the Landau functional with both experimental and ab-initio measurements as described below. The 
continuity of the first and second derivatives of the Landau functional are preserved which is sufficient for the present 
case. 

The next three terms with ^2,^3,^4 amplitudes are the contribution of the star (1 0). The four last terms 
associated with B2, -B4, B\, Bq are the contribution of the concentration waves {Qi}i. The expansion to the sixth 
order is required to obtain the linear stability of both DO22 and LI2 phases for the same range of concentration. 
The {-K4} and {B' 4 } coefficients couple concentration waves with different orientations. They control the amplitude 
of the potential barrier between the minima that correspond to ordered phases. The K3 term is the amplitude of the 
coupling between the waves that belong to the same orientational DO22 variant, (Ki,Qi) and equivalents. 

In fact, the precise form of the Landau functional term by term has not a direct influence on the mesoscopic 
microstructure providing that the functional is globally invariant with respect to the space group of the fee lattice. 
The important ingredients are the excess of free-energy associated with interfaces and long-range elastic interactions 
between domains. Therefore, we keep the lowest order coupling term between the LRO parameters associated with 
(Ki, Qi) , i. e. , terms of the form 7^77,:. 

The energy scale is fixed by the -Fb coefficient. The parameters of the adimensional free-energy density /l(co, 77, 7) 
are adjusted to fit the required qualitative thermodynamical properties. In Fig.l is plotted versus concentration Cq 
the free energy density Fq.Jl minimized with respect to the LRO parameters. The three types of minima correspond 
to the disordered phase, the DO23 (or DO22) and LI2 structures. The common tangent constructions determine the 
regions where ordered phases may coexist with the solid solution. The concentrations CD0 23 ,a and CD0 23 ,b are the limit 
of the region where coexist both disordered solid solution with the DO23 phase. Similar quantities can be determine 
for the coexistence of solid solution with LI2 phase: ci,\ 2 , a and C£i 2 .f,. The concentration CLi 2 . a and Coo 23 ,a are the 
solubility limit of the LI2 and DO23 phases, respectively. 

In a non-uniform system, the order parameters 77^ and 7, and the local concentration cq are spatially dependant. 
Within a phase-field approach, these parameters vary continuously thought the system. The energy excess due to the 
interfaces of precipitate, is expressed as a continuous Hamiltonian of the parameter fields and their first derivatives. 
To the lowest order, this leads to the following Ginzburg-Landau (GL) free energy density: 

Igl = F .[f L (c , V n) + {A c |Vc | 2 + ]T \ v V 2 m + A 7 V 2 7l }] (2.4) 

i 

where the A coefficients are the weight of the gradient terms. Total energy is given by Fql — J fGhdV. For the 
numerical implementation, we introduce a discrete space which is a cubic sub-lattice with unit cell of linear size 
d. This length must be large enough to justify our continuous approach and define the scale of one pixel in our 
simulations. The total free energy can be expressed as a discrete sum Fgl — d 3 Jcl where L represents the set 
of the sub-lattice nodes. 

We now describe the physical requirement we use to adjust the Ginzburg-Landau functional. First, cli 2 o, and 
CD0 23 a (see Fig.l) are the solubility limit of both LI2 and DO23 if each phase is supposed to coexist alone with 
solid solution. The cuo 22 a is given by the uncoherent phase diag rEim Enid CL\ 2a has been estimated by the measure 
of the lattice parameter of the solid solution by X ray using Vegard's law ||. For Zr at temperature T — 425°C, 
c D o 22 a = 0.0308 and c L i 2 a = 0.0426 at.% §. 

Second, to analyze the competition between the metastable LI2 and stable DO23 structures during the precipitation 
and aging processes, the free energy difference between both phases is an important quantity. In order to estimate 
this difference we refer to the theoretical studies of the formation energies obtained with ab-initio electronic structure 
methods j| In Q, it is found that ground-state of Al 3 Zr is indeed -DO23, and its structure is stabilized by the 
relaxation of the atomic positions inside the elementary motif. The difference of energy between LI2 and DO23 is 
found to be 5li 2 — 0.863 x 10 8 J/m 3 . In term of a phase field approach, this corresponds to an uncoherent energy 
at zero temperature. We assume that for low enough temperature, the free energy difference Sli 2 do not vary 
strongly. Therefore, the ab-initio quantity Sli 2 is used to fix the scale Fq of the free energy (Eq. |2.4D through the 
relation Fq = 8li 2 /5li 2 where 5li 2 is the energy difference between the corresponding LI2 and DO23 minima of the 
adimentional free energy which has been minimized taking into account the conservation of the Co concentration 
parameter that yields the common tangent construction (see Fig.l). 

Finally, another very important feature of the GL functional is the excess of energy of the interface between 
solid solution and precipitates, noted I with o = {LI2, D02s}- These quantities play a role in nucleation and 
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growth process. The interface energy of ordered precipitates in the solid solution can be measured experimentally 
at very low supersaturation. The Lifshitz-Slyozov- Wagner theory []lo|] gives the interface energy as a function of the 
diffusion coefficient which is physically measured. Unfortunately the few we found in literature (see &) about such 
a measurement is not satisfactory as the interface energy measured for LI2 precipitates Ili 2 is hundred times larger 
than the usual values. Thus we choose to estimate Ili 2 to a value of 10mJ/m 2 which is the order of magnitude 
of interface energies measured in aluminum compounds. As we did not find measurement in literature concerning 
interface energy of DO23 precipitates Ido 23 , we choose for Ido 23 a similar value to that of Ili 2 because there is no 
physical reason these two quantities to differ by an order of magnitude. We adjust the interface energies of the ordered 
phases {0} of the adimensional GL functional I such as F dI o = I . It implies Ido 23 /Ili 2 — Ido 23 /Ili 2 - With some 
difficulties, we managed to adjust the GL functional such that both Ido 23 and Ili 2 have the same order of magnitude. 
The Fq being fixed by the second criteriion stated previously, it imposes d — I Q j ' I Slx 2 /$li 2 so we can define the 
scale of our simulation. In order to investigate the nanometer scale we choose d = lnm, then the / functional must 
verify I /Sli 2 ~ 0.26. The only way we found to satisfy the previous criteria is to introduce the non-symmetric 
power n± for the term (c — c\) n± in the GL functional. Finally we obtain cuo 22 a = 0.0308%, cli 2 o, — 0.5% and 
Fo.d.lLi 2 = 8mJ/m 2 , Fo.d.Iz)o 22 = 9mJ/m 2 . The pixel of our simulation represents a cube of size d = 0.5nm. 



III. MICROELASTICITY CONTRIBUTION TO FREE ENERGY 

As described in |]lT|| , the elastic energy E e i is calculated assumsing that the local strain (e^) induces a relaxation 
that is calculated by setting a small volume dV of the bulk to the mechanical equilibrium. It is supposed that the 
time needed to reach the mechanical equilibrium is negligible compared to the typical diffusive time of the ordering 
process. The key point of the phase field theory for alloys is that the stress free strain tensor can be expressed as a 
function of the local LRO parameters and local concentration cq(R). 

The geometrical operation to transform a cubic unit cell of the solid solution with lattice parameter a into the cubic 
unit cell of the LI2 phase with lattice parameter o.l\ 2 is given by the tensor: 

= s kl (aLi 2 - a) /a (3.1) 

where Ski is the unity tensor. The lattice parameter a^i 2 has been measured for a perfect LI2 phase j^] i. e. , with 
stoichiometry 0.25 at. % Zr and a is extrapolated from the lattice parameter ao of pure aluminum Al, using Vegard's 
law. 

We note auo 23 and boo 23 the lattice parameters of the tetragonal DO23 phase at stoichiometry 0.25 at. % Zr. 
The geometrical operation to transform a cubic unit cell of the solid solution into the tetragonal elementary cell of 
the DO23 phase is given by the tensor e^° 23 (p) if the orientational variant corresponds to the association of the 
(K p , Q p ) waves with p = 1,2 or 3: the cell is dilatated in either (100), (010) or (001) direction. Here we choose as an 
example the orientational variant associated to (K\,Q\) and the cubic unit cell is dilated in the direction (100) so as 
i>oo 23 >a> a,D0 23 - 

_- 1 1 

e DO ™(l) = gggg - [ t I (3.2) 

( 

where t = {ano 23 — &)/(&do 2 3 — a). Tensors 23 (2) and 23 (3) derive from e^° 23 (l) by permutation of the 
diagonal coefficients. For the general expression of the local strain, we propose the following form (see 

3 

e° kl (R) = e°?(0) MR) + E e ™ (P) MR) (3-3) 

P =i 

where ijj p (R) = (7 p (i?)) 2 , and ipo(R) = (cq(R) —c). The tensor coefficients e^(p) are chosen in a such way that 
e ki(R) ~ e ki 2 if a re l axe< i LI2 inclusion is at the position R and ^(R) — e^ 023 ^ if a DO23 inclusion is in the same 
position with the orientational variant p. The Eq.3.3 is rewritten in a compact form: e°(R) = X)p=o e °°(p) ''Pp(R) 
where the p indice varies now from to 3. The strain- induced elastic energy can be computed following the model 
presented in ]ll[ which gives: 

Ed = \zZ I [B Pg Wp^ q dK 3 (3.4) 
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where ip p is the Fourier transform of the function ip p . There 

B p . q (K) = Xijki eg?(p) e°?(g) - b pq (K) 

and b pq (K) = of°(p) ^ Gj^ if; ^(ff) where Ayfc/ is the elastic tensor and af?(p) = Kjki^ip)- The tensor Gjk is 
the elastic Green function. The tensor Xijki is assumed to be homogeneous in space and the simulations are realized 
with elastic coefficients of aluminum [jL2| . 



IV. KINETICS OF THE PHASE TRANSFORMATION 



The total energy is given by the sum F = Fql + E e i (Eqs.( |2.4|) and (3.4)). At mesoscopic scale, the kinetics of 
the phase transition is well described by a phase field method (see |l3|]). In the context of a phase field approach, 
the local composition Co is a conservative order parameter and thus its time evolution is driven by the Cahn-Hilliard 
equation: 

and for the non-conservative LRO parameters the kinetic is given by 

f = -L,g +,,*(*,*) (4.2) 



9 -i = -L^ + r , n.n 



where v's are stochastic terms. To simulate thermal fluctuations, it is useful to introduce the Langevin noise which 
consists in assuming a white space-time noise for the stochastic terms and no cross-correlation between each other. 
Numerically, the r andom fu nctions v's are implemented with a gaussian probability density 

The set of Eqs.[(4.1)-(4.3)] is the so called Time-Dependant Ginzburg-Landau equation ||| and they can be derived 



from the microscopic Onsager equation with respect to the occupation probability of the solute atoms Q . Numerically, 
for simplicity, we choose L„ = L„ = L c = 1. So in order to deduce the approximative time unit of our simulations, 
our numerical time must be divided by the numerical value of the diffusion coefficient. The Langevin noise provides 
a very primitive description of the thermal fluctuations. All though it is the simplest and the less controversal way to 
implement the temperature in the simulation. 



Equations [(4.1 )-(4.3)] are integrated on a system that is invariant along z-axis to save computation time. The 
initial state of the system is a uniform solid solution at concentration Cti 2ia < c = c < 0.25 and all the set of LRO 
parameters are set to zero value, which represents an unstable disordered phase. The initial time of our simulation 
can be considered as an instantaneous quenching of the material. 

To represent the microstructures we choose to color the DO23 precipitates with either blue or red depending on 
their translational variant. The four translational variants of LI2 are allocated to four other paler colors. The 
disordered phase is colored with black. The gray scale is sufficient to distinguish the different types of precipitates 
if the translational variants are forgotten. On Figs. [(??)-(??)], we present the dynamics of the phase transition from 
disordered solid solution to a microstructure with ordered precipitates embedded into a disordered matrix. The 
different sequences are realized for different average compositions c and different temperatures T . 
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First, one notes the specific rod shape of the DO23 precipitates. The tetragonal symmetry combined with a large 
misfit (&DO23 — I® which involves a large intrinsic stress is well known to induce such a pattern [fi"l| . On the very 
last picture of Fig. 2 the facets of the DO23 inclusions correspond to the habit planes with orientation of around 20 
degrees with respect to the (1 0) directions. For any couple of external variables, namely the temperature T and 
the composition c the late stage of the kinetics is a microstructure which contains exclusively the DO23 precipitates 
embedded in a disordered matrix. Nevertheless LI2 structures may appear in the early regime of the dynamics (see 
Fig. 3). As the LI2 ordered precipitates involve a weak misfit (q.li 2 — a)/a compared to (&DO23 — <^)A*i their shape is 
spherical. These spherical inclusions nucleate only at low temperature T and low Zr concentration c. At low enough 
values of T and for saturation c close enough from the solubility limit, the microstructure contains exclusively the LI2 
precipitates inside the solid solution. Since the Ll 2 phase is metastable, any grain of this phase should not resist to 
thermal fluctuations and thus no LI2 precipitates should grow in the microstructure. It is actually what is observed 
if the elastic energy E e i is neglected in our simulations. However, in the limit of low Zr saturation, the kinetics 
of the transition drives the system to a transient LI2 microstructure which is favored by its intrinsic stress which 
is weaker than that of DC>23- If the temperature T is not too high, i. e., T < Tl\ 2 m 500K (see Fig. 2), the LI 2 
precipitates can even grow by consuming the very few solute atoms contained in the disordered matrix. This implies 
that the solid solution becomes poor in solute and therefore the nucleation of DO23 precipitates is no longer possible 
in the disoredered matrix. In fact, if temperature is lower than Ti,i a , the LI2 microstructure may survive for a time 
much longer than the computation time. Nevertheless, a gradual increase of the temperature reveals the process of 
nucleation of the DO23 precipitates (see Fig. 2). A remarkable result is that nucleation of the stable phase occurs at 
the interface of the LI2 precipitates with the solid solution where the local concentration cq is high enough. One 
note that the preferential growth of DO23 precipitates occurs at temperature higher than T^i 2 that demonstrates the 
robusness of the metastable Ll 2 microstructure with respect to thermal fluctuations. 

Similar simulations at larger saturation c show that both the DO23 and LI2 precipitates nucleate and may coexist 
in the microstructure (see Fig. 4). In that range of c, the DO23 inclusions can nucleate at places different from LI2 
precipitates because the solute atoms have not been consumed by the growth of LI2 inclusions. Once saturation 
of the solid solution is locally dried up, the DO23 precipitates grow to the expense of the LI2 grains via the solute 
diffusion throught the matrix. Then the persistent LI2 precipitates are localized relatively far from the DO23 grain. 
For a given temperature, if the average concentration c is increased again, the grains with different phases nucleate 
in neighboring regions of the supersaturated solid solution and the DO23 inclusions absorbe the Zr matter of LI2 
inclusions. 

One remark on Figs.[(2)-(4)] that the orientational variant of DO23 which is combination of waves {K^^Q^) is 
inhibited. It is because the precipitates with such variant cannot relax their elastic energy because of z-invariance. 
Qualitatively it is not a problem as there are still two orientational variants for the DO23 phase. 

V. CONCLUSION AND PERSPECTIVE 

The present paper treats of the specific case of the interplay between the LI2 metastable phase and the DO23 
ordered ground-state during the order-disorder transition in Al 3 Zr alloy. It is proved that for a sufficiently low 
temperature and weak solute saturation the metastable phase nucleates before the stable phase. It is the result 
of the dynamics of the phase transition that is deeply influenced by the microelasticity induced by the strain of 
precipitates. The metastable ordered phase with a weak internal stress may be favored with respect to the ordered 
ground-state which induces a much larger strain. The strain-induced elasticity may play a role at the early regime of 
the dynamics of the phase transition. Furthermore, depending on external variables, temperature and composition, we 
found different kinetics for the vanishing of the LI2 precipitates. At low temperature and low saturation of solute, the 
DO23 precipitates grow preferentially at the interfaces of LI2 inclusions with solid solution. At higher temperature or 
equivalently higher Zr saturation, the DO23 precipitates nucleate into the solid solution. Once DO23 precipitates have 
nucleated, they grow at expense of the LI2 structure that disappears. Actually, such phenomena have not yet been 
observed experimentally. With this respect, the phase field method can be considered as a predictive method, though 
experimental confirmation is now required. On that trail, microscope analysis are programmed in the Laboratoire 
d'Etude des Microstructures (ONERA). 

It is well known experimentally that the LI2 microstructure has better mechanical properties than DO23. As 
the LI2 is metastable, the degradation of the mechanical properties with increasing temperature cannot be avoid. 
Nevertheless our results allow to hope that it is possible to increase the robusness of the LI2 microstructure playing 
with elastic interaction. We expect our study will contribute to the understanding of precipitates formation and to 
improve the control of the alloys synthesis. To that purpose, the phase field method we use, can be extended for other 
alloys with similar phase transition as T13AI or Pd^V . 
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Finally our calculations are only valid for coherent sample at the nanometer scale. It is possible to investigate 
non-coherent effect introducing dislocations in the phase-field method as it is discribed in both E3] jl6| . It is of great 
interest to perform such simulation in the case of Al^Zr where discontinuous precipitations and dislocations modify 
strongly the precipitation process J^,[^| and play an important role in macroscopic properties of materials. 
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